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We study the contribution of massive douiinantly sterile neutrinos, A*', to the Lepton Number 
and Lepton Flavor Violating (LNV, LFV) semileptonic decays of r and _B, D, .ft'-mesons. We fo- 
cus on special domains of sterile neutrino masses nijv where it is close to its mass-shell. This 
leads to an enormous resonant enhancement of the decay rates of these processes. This allows 
us to derive stringent limits on the sterile neutrino mass mjv and its mixing matrix elements 
(**^ |l7ev|, |t/^v|, IC^evf/rivl, It/rvt/rivl, by Carrying out a joint analysis of the existing experimental 

y—{ bounds on the decay rates of the studied processes. In contrast to other approaches in the litera- 

(^ ture, our extraction of these limits is free from ad hoc assumptions on the relative size of the sterile 

C\| neutrino mixing parameters. Special attention was paid to the limits on meson decays with e^e* 

K^ in final state, derived from non-observation of Oz^/3/3-decay. We point out that observation of these 

^ decays may, in particular, shed light on CP violation in the neutrino sector. 

^ PACS numbers: 12.15.-y,12.60.-i,13.15.+g,13.20.-v,14.60.Pq 

- Keywords: bottom mesons, neutrino, decay rates 



Oh 



> 

O 



in 

o 
o 



X 



I. INTRODUCTION 

Neutrino oscillation experiments have shown conclusively that neutrinos are massive, although very light, particles 
mixing with each other. Moreover, neutrino oscillations is the first and so far the only observed phenomenon of 
lepton flavor violation (LFV). Theoretically LFV is only possible if neutrinos are massive. Being electrically neutral 
neutrinos can be either Dirac or Majorana type particles. Majorana neutrino masses violate total lepton number 
conservation and, therefore, can induce lepton number violating (LNV) processes. Up to date such processes have 
not yet been observed experimentally. Searching for LNV processes is a challenging quest, pointed to probe the 
nature of neutrinos and answer the question about wether they are Majorana or Dirac particles. However it is well 
known that in Standard Model (SM) extensions with very light and very heavy neutrinos, the rates of both LFV and 
LNV processes are so small that their experimental observation turns out to be unrealistic. Perhaps the only lucky 
exception is neutrinoless double beta decay (Oz^/3/3). The experiments searching for this LNV process are believed to 
reach even such small neutrino masses [B [2] as those which are relevant for neutrino oscillations, and to probe the 
heavy Majorana neutrino sector up to the masses of several TeV 0|lj. 

On the other hand there are other extensions of the SM in which both LFV and LNV effects could be significant for 
various reasons. In particular there may contain new massive particles like moderately heavy Majorana neutrinos, as 
in the low-scale versions of the see-saw scenario, or other new particles such as neutralinos, heavy Majorana particles 
mixing with neutrinos in the framework of SUSY models with R-parity violation. These new heavy particles may 
lead to observable effects beyond the light neutrino sector, since they could manifest themselves indirectly via their 
virtual contribution to processes involving ordinary particles. If these new particles are not very heavy lying within 
the mass reach of the running of forthcoming experiments, they can be searched for directly among the products of 
colliding particles. Also they may contribute to some decays as virtual, but nearly on-mass-shell states, leading to 
huge resonant enhancement of the corresponding decay rates. 

Here we consider an extended see-saw scenario, including n species of SM singlet right-handed neutrinos z/^ ■ = 
(j/^;^, ...r/^jj), besides the three left-handed weak doublet neutrinos iy[^ — {vl^^vl^.v^^). The general mass term for 
this set of fields can be written as 

i^AlM.- + h.c. = -Ic^',^)!"-^^ f/ Vn+h-c- (1) 

2 2^ "^^ '''\M% Mr ) \v'^ )^ ^' 

-.3 n 

= -2C^m^,'i^i^, + J2^^j'i^''j)+^-(^- (2) 

Here Ml, -Mr are 3x3 and nxn symmetric Majorana mass matrices, and Md is a 3 x 7i Dirac type matrix. Rotating 
the neutrino mass matrix to the diagonal form by a unitary transformation 

[/r7W('')C/ = Diag{m,i,---,m,3^„} (3) 

we end up with 3 + n Majorana neutrinos with masses m^j , • • • , fnv^^,^ ■ In special cases there may appear among 
them pairs with masses degenerate in absolute values. Each of these pairs can be collected into a Dirac neutrino field. 
This situation corresponds to conservation of certain lepton numbers assigned to these Dirac fields. 

Neutrino interactions are represented by the SM Charged (CC) and Neutral Current (NC) Lagrangian terms. In 
the mass eigenstate basis they read 

-^ = 4 E U^^ '^^'P^^^ ^i^ + ^^fe^ E U.,K, i^^i'Pli^, Z„ (4) 

where I — e, fi,T and i — 1, ..., n -I- 3. 

In this scenario the light neutrinos have the mass scale M^/Mr, and as a result there should be also heavy 
Majorana neutrinos (N) with mass scale Mr. In particular, they may include moderately heavy neutrinos N in the 
MeV-GeV domain, which we are going to consider in the next sections. If they exist they may contribute to some 
LNV and LFV processes as intermediate on-mass-shell states. This would result in enormous resonant enhancement of 
their contribution to these processes and would allow one to set stringent limits on heavy neutrino masses and mixing 
from non-observation of the corresponding processes. The presence or absence of these neutrino states, which are 
conventionally called sterile neutrinos, is a subject to be revealed by experimental searches. Having these arguments in 
mind we study neutrino states N with masses rriN in this region of values and examine their phenomenological impact 
on the LFV and LNV decays of t and some mesons. The phenomenology of sterile neutrinos have been studied in the 
literature in different contexts and from different points of view. Their resonant contributions to r and meson decays 
have been studied in Refs. [S], [B], [7]. Another potential process to look for sterile Majorana neutrinos is like-sign 



dilepton production in hadron collisions [5], [5], [10] i [H]- Possible implications of sterile neutrinos have been also 
studied in LFV muonium decay and high-energy muon-electron scattering [12, . (for a recent review of sterile neutrino 
phenomenology see [13]). The present paper contributes to some still uncovered aspects of the phenomenology of 
sterile neutrinos. 

The paper is organized as follows. In the next section In] we present decay rate formulas for t and meson LFV and 
LNV decays and their reduced forms, in the resonant domains of neutrino masses. In sec |III| we discuss theoretical 
uncertainties in the calculation of the total decay width of sterile neutrinos and present our "inclusive" approach based 
on Bloom-Gilman duality. This approach allows one to avoid large uncertainties related to the heavy meson decay 
constants, which are typical in the conventional channel-by-channel approach. In sec IV we discuss the extraction of 
the limits on the active-sterile neutrino mixing matrix elements |t/eAr|, |t^/jjv|, WenUtN , {U^^nUtnI and present the 
corresponding exclusion plots. In sec. IV] we discuss the Oi'/S/J-decay constraints on sterile neutrino parameters and 
their possible implications for CP-violation in the neutrino sector. Here we also give predictions for the rates of some 
LNV and LFV decays of K, D, Ds, B, i3c-mesons. Sec. VI summarizes our main results. 



II. DECAY RATES 

In what follows we analyze the above mentioned resonant contribution of heavy sterile neutrino to the semileptonic 
LNV and LFV decays of t and the pseudoscalar mesons M — K,D,B: 



/=F7r±7r-, 



M+ ^ TT^lfli 



(5) 



Lowest order diagrams for the case of meson decays are shown in Fig. [T] There is only one tree-level diagram with an 
intermediate Majorana or Dirac neutrino, shown in Fig. RJa), contributing to LFV decays. For LNV decays mediated 
by Majorana neutrinos, in addition to the tree-level diagram Fig. [Ha), there appears a two- loop diagram shown in 
Fig. [ijb). The tree- level diagram is known fT to dominate in these processes ((sl). As will be seen, in the studied 
domain of the heavy neutrino mass niN, the two-loop diagram in Fig. WCb) is absolutely negligible. An important 
point is that the calculation of the tree-level diagram. Fig. ^a), does not require knowledge of the hadronic structure 
needed for the diagram in Fig. [lib). 
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FIG. 1: The lowest order diagrams contributing to the semileptonic meson decays. 
The decay rates of the studied processes are given by the expressions: 
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LFV : T{M+ -^ tt^I'IJ) = 
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Here the unitary mixing matrix Uij relates v[ = Uiji^j, the weak v' and mass ly neutrino eigenstates. In Eq. (rol) the 
factor Kij = 1 — Sij/2 takes into account a combinatorial factor 1/2 for identical final leptons. In the above equations 
we introduced the following functions: 



Hm{zi,Z2) 



c'"^ [hUz) - xlh,+{z)] [hMz) - hU^)] 



(8) 
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— I hi^{zi)hj^{z2) + xl [h.,^{zi)hj^{z2) + xj + xf\ - 

M 

[hi^{zi) + /ij-(z2)] [hi-{zi)hj+{z2) + hi+{zi)hj^{z2)] 



with 



c^^3{z) = Ai/2(z, xlxl)\^/\z, xl 1), h,±{z) = z±x 



„M 



G\ 
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nfMK,\yudV\VM\\ Xix,y,z) ^x'+y' + z'- 2xy - 2yz - 2xz 



(9) 
(10) 



and the dimensionless variables Xi 



mijm^ 



TOtt/TOjj^, where m/^ and m^^j are the masses of the charged 



leptons li = e,fj,,T and the initial meson M respectively. Numerical values of meson masses ttim, decay constants /j^, 
and the CKM factors V^/ are specified in Table I. The integration limits in Eqs. ^, (l7| are 



m\{x^ + xif 



<(i-^0', 



(11) 
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with y = sxjrr? . 

If one assumes that neutrinos are separated into light v^^ and heavy N\^ states, with masses rriyi « 
<< m]^i_, then the branching ratios '5] turn out to be extremely small: 
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3 light neutrino scenario. 
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3 light + 1 heavy neutrino scenario 
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These values are far beyond experimental reach. 

On the other hand, if we assume that there exists at least one massive neutrino N with mass m^ in the range 



< rriN << niN < max 



(14) 



the situation changes drastically. Then the s-channel M-decay diagram in Fig [TJa) blows up because the integrands 
under the single integrals in Eqs. ([6]), ([7| have a non-integrable singularity at s = m% corresponding to an on-mass- 
shell intermediate neutrino. Therefore, in this resonant domain the total decay width of heavy neutrino r^r has to be 
taken into account. This can be done by the substitution m^ — >■ rnj^ — {i/2)Tiq. As will be seen in sec. IV the heavy 
neutrino width Fjv is very small in the resonant domain, F^v ~ 10~^ MeV, [SlfB]. Therefore, the neutrino propagator 
in the single integrals of Eqs. (|6|, (|7l has a very sharp maximum at s = raj^. The double integrals, being finite in 
the limit T^h = 0, can be neglected in the considered case. Thus, with good precision we obtain from Eqs. ([6]), (l7| the 
following decay rate formulas for a nearly on-mass-shell intermediate neutrino: 



LNV : F""(A/^ 
LFV : r'^"(M^ 



^ K^) ~ KijTT{Glj{zM) + G\,i{zm)) 



,mN\Um\^\UjN\^^ 
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(15) 
(16) 



with zm = ("T-Af/'TiAf )^- Similar arguments apply to tau decays r —> Z+tt 7r+ "^j which in the neutrino mass domain 
m,r + '71; < JTijv < "Ti^ — JTiTr lead also to resonantly enhanced decay rates: 



where 
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(G^^/32)(7r)-V>^|K. 



(18) 



with z^ = (rnj<i jra-r)'^ . We will use the above decay rate formulas Eqs. ( 15 1-( 18 ) in our analysis of the sterile neutrino 
contribution to t and meson decays. 



III. NEUTRINO DECAY WIDTH 

Heavy sterile neutrinos A^, being mixed with active neutrino flavors, can decay in various final states depending on 
their mass ni]^. These decays are represented by purely leptonic TV — >■ l\l2V., ?>v and semileptonic N — >■ IH modes, 
where H — M, _B, .. are hadronic states represented by mesons and baryons. These decays proceed via CC and NC 
interactions of the SM given by the Lagrangian Q. The total decay rate. Fat, of heavy neutrinos is conventionally 
calculated in the literature in the channel-by- channel approach [51 1^, in which one sums up the partial decay rates 
of all the leptonic and two-body semileptonic decay channels open for a given value of ttiat. In the present paper 
we apply another approach. This is the inclusive approach we proposed in Ref. [5], in which we approximate the 
semileptonic decays of the heavy neutrino N by its decays into quark-antiquark pairs N -^ l{i')qiq2, as suggested 
by Bloom-Gilman duality |14j . This implies that in an average over sufficiently wide range of the invariant mass of 
the final hadronic state H, the sum of all the open decay channels N — > l{v)H is approximately equal to the rate of 
N — )■ l{v)qiq2- In comparison with the former approach the latter does not require information on the parameters of 
the final mesons, such as masses and decay constants, some of which are poorly known for the meson states starting 
from 77' (985) and heavier. The inclusive approach is supposed, according to the duality arguments, to take into account 
all the semileptonic channels and therefore, in this case Fa? should be larger than in the channel-by-channel approach, 
in which some hadronic states are neglected. We apply a simplified version of the inclusive approach, neglecting 
perturbative and nonperturbative QCD corrections to the tree-level quark production diagram. This leading-order 
approximation is expected to be reasonable for miq > > A sa 200 MeV. At lower masses a more viable approach would 
be to relate by dispersion relations the semileptonic N decay rate to the imaginary parts of the W and Z self-energies 
n(s), in analogy to the approach applied in the literature for the r — > v+hadrons inclusive decay [TS]. However for our 
rough estimations we do not need this more sophisticated treatment and will use the above mentioned leading-order 
approximation. 

Here we summarize the partial decay rates for the inclusive approach [B] , including leptonic and semileptonic decay 
modes of the heavy sterile neutrino N. In the latter case the final hadronic states for low neutrino masses m^ is 
represented by the lightest mesons while for larger uin by gg-pairs. The list of decay rates is as follows: 
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Table I: Masses mp^v and decay constants /p,v of pseudoscalar P and vector V mesons. 
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775.8 


220 [13] 


Kd 


K^ 


493.7 
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222.6 |24j 
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jj,± 
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Dt 
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B± 


6277 


399 m] 


Kt 


LU 
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tt" 
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r(iV^/r^+) = \Ui,N?^^llv\Vv?Fv{yi,.yv)^\Ui,N?T^'''\ 
T{N^vi,qq) ^ \Ui,N?—^m%Fz{mM)[g^L9Rh{y,.,^,yq,yq) + 

+ {{g'Lf + {9Rr))h{y.,,,yq,yq)] ^ |f//,^pr(''«^^ 



(24) 
(25) 
(26) 
(27) 



Here we denoted yi — rrii/mi^ with mi = mi,mp^mv ,mq. For the quark masses we use the values 
rriu « md = 3.5 MeV, m^ = 105 MeV, m^ = 1.27 GeV, m{, = 4.2 GeV. In Eqs. ^, (IStI) we denoted u = m, c, i; 



d = d, s, fe and q = u,d, c, s, b, t. The SM neutral current couplings of leptons and quarks are 

gi = -l/2 + sin2 0^, gl^l/2~{2/i)siu'0w, g^ = -1/2 + (1/3) sin^ ^vy, 
5^ = sin^ 0^^, g]^ = -(2/3) sin^ 6w, 9r, = (1/3) sin^ 6^^ 

The corresponding NC couplings of mesons are given by 

Kv — sin^ 0w/3 for p^,(jJ, 

Kv = -1/4 + sin2 6*^/3 for K*°,K*°,(l3, 

Ky = l/4-2sm'^ew/3 for D*°,D*'^, J/iP 



(28) 



(29) 



The kinematical functions are 



/i(x,y,z) = 12 / ^{s~x'-y'){l + z'-s)X'/'{s,x^y')X'/^l,s,z'), 



{l-xf 



l2ix,y,z) = 24yz 



ds 



il + x'~s)X'/'is,y',z')X'^\l,s,x'), 



{y+zy 



Fp{x,y)^X'/'il,x\y^)[{l + x'){l + x'-y^)-4x^] 
Fv{x,y) = Al/2(l,x^y2)[(l _ ^2)2 ^ (i + ^2y _ ^y4 



(30) 



The total decay rate T^ of the heavy neutrino N is equal to the sum of the partial decay rates in Eqs. (19 1- (27 1, 
which we write in the form: 

Tjv = ^ [mnN^iYH+) + mnN^i^itiyi,)+ (3i) 

where we denoted the hadronic states H^ = P^,V^,du,su,dc,sc and H^ = P'^,V'^,qq. We introduced the factor 
r^AT = 2 for Majorana and 77JV = 1 for Dirac neutrino N. Its value lyA? = 2 is related with the fact that for Majorana 
neutrinos both charge conjugate final states are allowed: and N — >■ Zj"Z^z/;2, ltl2'^h ^^^ -^ ^ PH^. For convenience 



we write Eq. (31) in the form: 

Tpf = ae{mN) ■ \UeN\^ + a-iiimN) ■ {U^j.nI'^ + ar{mN) ■ \UrN\^ (32) 



where 



ai (m^) = r(™) + r(3-) + ^ [r^hhu) ^ ,^^r("^^)) , (33) 



and I, I2 — e, ^, r. In the inclusive approach the hadronic contribution is calculated as 

r(™) ^ 0(^0 - mjv) E (r^'^'^^ + r^''^^ + m^^'^^ + ^r^'^^) + 0(m^ - ^o) E ('y^r^''"^) + r^'"'^)) (34) 

P,V u,d,q 

The parameter /ip denotes the mass threshold from which we start taking into account hadronic contributions via qq 



production. In our numerical study we use the mass /xq = rriK = 493.7 MeV. Thus in the first term of Eq. (34) we 
take into account only P — tt^ , tt*^. At this rather low threshold (/io) the QCD corrections become significant, but it 
is reasonable to expect that the tree-level contribution still dominates the semileptonic decay rates since /iq > ^qcd- 
We believe that the accuracy of the above presented inclusive approach is sufhcient for estimations of limits on the 
parameters of sterile neutrino in this range of masses. In this respect it is worth remembering that in the channel-by- 
channel approach uncertainties related to the heavy meson decay constants Jm and other hadronic parameters are 
large and theoretically not well controllable. 

IV. LIMITS ON STERILE NEUTRINO MASS AND MIXING 

In the literature there are various limits on the mass rriN and mixing UaN (with a — e, fj,, t) of a sterile neutrino 
N, extracted from direct and indirect experimental searches [TB] for this particle, in a wide region of its mass. A 
recent summary of these limits, extracted from the corresponding experimental data, is given in Ref. |13j . In Figs. 
m3] we show the exclusion plots from Ref. |13j together with our exclusion curves, derived from a joint analysis of 
semileptonic LNV and LEV decays of K, B, D, D^-mesons and r. 

As seen from Eqs. (15), (17) and (32), the decay rates of these processes depend on all the three Ue^f^^rN mixing 



matrix elements. In the literature it is a common practice to adopt some ad hoc assumptions on their relative 
sign in order to extract limits on them from the experimental bounds on the corresponding decay rates. These 
assumptions reduce the reliability of the obtained limits. In what follows we carry out a joint analysis of the above 
mentioned processes, without any additional assumptions of this sort. We apply a numerical Monte Carlo sampling 
in a parametric space 

\UeN\,\Uf,N\AUrNlmN, with 0<|C/iAr|<l, mNeRes{M,T) (35) 

taking into account the existing experimental bounds on the branching ratios of LNV and LEV decays: M —>■ lln, 
T — > /tttt. Here Res{M,T) denotes resonant regions of sterile neutrino mass corresponding to the studied meson 
M = K, B, D, Dg and r decays. The experimental limits on their branching ratios were taken from Ref. [16 , and are 
shown in Table II. The obtained exclusion curves for \UeN\ and \U^n\ are shown in Figs. l2J [sl together with other 
existing limits. Figs. Hshow the exclusion curves for |C/eAfC^rAr| and |[/^Arf/,- wi- 
lt is sometimes useful to have analytic expressions for the limits on sterile neutrino mixing. Such expressions for 
the case of |C/eArP and |t/^Arp, without ad hoc assumption on the relative size of different matrix elements |C/iArp, can 



be derived directly from Eqs. (15l - (17) and (32). Using the fact that |t/^Ar| < 1, |f/rAr| < 1 one finds 



WeM" - \U,N?F,,[M)a, - F,e{M){a^ + a,) < (36) 

\U^n\^ - |t/^jvpi^Mp(M)a^ - F^^iM)ia„ + a,) < 



Solving these inequalities one gets limits 



\UeN\ < Z 



\U^n\^ < 



F^^{M)a^ 




Fee{M)ae 



-^ ee \^fj 



F^.^.{M)a, 



+ F^^{ae + ar) 



valid in the resonant regions of decays M — >■ eeTr and M — > /i/i7r. We introduced the notation 



F^e(M) 






Fu{M) 



cMTmiNG'HzQ)' 



Fi{r) 



c"iTmNG\zr) + 6G"='"(zo))' 



(37) 
(38) 

(39) 



in agr eem e nt w ith Eqs. (15)-(18|. Here the T'^^p are experimental upper bounds on the rate of the indicated decays. 
Eqs. (37), (38) allow setting limits on |C/eAfP and |C/^Arp, considering different LNV and LFV decays separately, and 



using the existing experimental bounds on their rates. 

Stronger limits can be extracted from a joint analysis of certain sets of LNV and LFV decays, if these decays have 
nontrivial intersection of their resonant regions in tti^t. Then for values of mjq in the intersection one can extract 
limits for Uin in an analytic form in the following way. Let us consider, for instance, the set of LNV and/or LFV 
decays: M — )• 7re^, M — )• Tree and r -^ inre. Then from Eqs. (15 1 - (17) and (31) we obtain 



ae\UeNr 

Fr{r) < 



\UeN\^ 
ra^\Uf_,N\^ 
\UeN\'\UrN\' 



- ar\Ur 



N\ 



Fet^iM) < 

ae\UeNr 



\UeNm^.N\' 
+ afi\U^N\^ +ar\UrN\ 



ae\UeN\^ +a^\Uf^N\'^ +ar\UrN\^ 

conibinig these Eqs. we find 

|C/eJv|' < aeFee{M) + a^F,^,{M) + arF,{T) 
Analogously, from the set M — > tt/x/x, A'I — )■ 7re/i and t — > tttt/i we get the limit 

iC/^ATp < aeF,^(Af) + a^F^^(Af) + arF^ir) 



(40) 
(41) 

(42) 
(43) 



Using in these formulas ( 42 ) , ( 43 1 , the experimental bounds on the branching rates of the corresponding decays [H] , 
one can analyze all the sets of decays of the above mentioned type for i?, Z?, Z?s, ^-mesons, and extract exclusion 
curves in |C/eArP ~ "f^N and |f/^Arp — m^ planes. The resulting exclusion curve for |L/eArp is limited to the mass range 
TOjv = 141 — 1637 MeV, while for |C/;^Arp to the range rriN = 246 — 1637 MeV. We do not present these curves since 
they are not very different in these limited mass regions from those derived in the Monte Carlo approach and shown 
in Figs. 2]|4 The same is true for Eqs. (37), based on the analysis of different decays separately. Thus the above 
presented analytical formulas (42), (43) or (37), considered in the corresponding regions of sterile neutrino mass ttiat, 



can be treated as a reasonable alternative to a more involving Monte Carlo analysis. 



V. LIMITS: COMPATIBILITY AND INTERPLAY 

Let us study some implications of the interplay and compatibility of the constraints on the same combinations of 
parameters derived from different processes. 



A. Impact of Oi^/3/3-decay limits. 

The exclusion curve from neutrinoless double beta decay {QujSp) in Fig. [2J despite the presence of some uncertainty 
in nuclear matrix elements (of a factor ^ 2), is so stringent that within this uncertainty it overrides all the other 
constraints except for the narrow region where the exclusion curve of the beam dump PS 191 experiment (TS] goes 
lower. However the PS 191 curve was obtained on the basis of various ad hoc assumptions, touching upon both data 
processing and their theoretical interpretation, making the corresponding limits not too firm to compete with the 
0;^/3/3-limits. In particular, the PS 191 limits can be evaded if one admits that the sterile neutrinos have dominant 
decay channels into invisible particles [131 [12] . 
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Table II: Experimental and theoretical upper bounds on the branching ratios of meson and t decays. The experimental 
bounds are taken from Ref. [16]. The theoretical upper bounds are shown in the square and curly brackets and are discussed 

in sec. V. 



Br{M -^ nhh) 


K 


D 


D, 


B 


Be 


M+ -^ 7r-e+e+ 


6.4 X 10"^" [1.4 X 10"**] 


3.6 X 10-^ [8.6 X 10"''] 


6.9 X 10-^ [1.0 X 10-^] 


1.6 X 10"'' [2.3 X 10""] 


{4.1 X 10-"} 


M+ -^ 7r-/i+^+ 


3 X 10"^ 


4.8 X 10-'^{3.1 X 10"''} 


2.9 X 10~^{3.7x lO"'^} 


1.4 X 10-''{3.4 X 10""} 


{5.4 X 10-*} 


M+ -^ 7r+//-e+ 


1.3 X 10"" [7.2 X 10-**] 


3.4 X 10"^ [2.5 X 10""] 


6.1 X 10"'' [3.1 X 10"®] 


1.7 X 10-^ [5.7 X 10-"] 


{1.0 X 10-"} 


M+ -^ ■K-ij+e+ 


5 X 10"^° [6.9 X 10"*] 


5 X 10"^ [1.7 X 10""] 


7.3 X 10"'' [2.2 X 10"''] 


1.3 X 10-*^ [4.6 X 10""] 


{ 8.3 X 10""} 


M+ -^ 7v+T-e+ 


- 


- 


{7.1 X 10-"} 


{8.4 X 10"^^} 


{1.4 X 10"''} 


M+ -^ Tv~r+e+ 


- 


- 


{7.1 X 10-"'} 


{8.4 X 10-12} 


{1.4 X 10"''} 


M+ ^ -K+T-fi+ 


- 


- 


- 


{2.5 X lO-'-*} 


{4.3 X lO-'^} 


M+ -^ n-r+fi+ 


- 


- 


- 


{2.0 X lO-''} 


{3.5 X 10-^} 


M+ -^ K-e+e+ 


- 


{4.0 X 10"^^} 


{5.3 X 10-"} 


{1.6 X 10""} 


{3.1 X 10-12} 


M+ -^ K+fj,-e+ 


- 


{9.8x10-^2} 


{1.3 X 10"} 


{4.2 X 10-"} 


{7.8 X 10-"} 


M+ -^ K-fi+e+ 


- 


{8.0 X 10"^^} 


{1.0 X 10-"} 


{3.4 X 10-"} 


{6.3 X 10-"} 


M+ ^ K-n+fi+ 


- 


{7.7 X 10-^} 


{1.0 X 10-^} 


{1.7x10-"} 


{2.8 X 10--'} 


M+ ^ K+T~e+ 


- 


- 


- 


{9.2 X 10-"} 


{1.6 X 10-"} 


M+ ^ K-T+e+ 


- 


- 


- 


{9.2 X 10-"} 


{1.6 X 10-"} 


M+ -^ K+T-fi+ 


- 


- 


- 


{1.4 X 10-"} 


{2.4 X 10-**} 


M+ -^ K-r+fi+ 


- 


- 


- 


{1.2 X 10-"} 


{2.1 X 10-*} 


Br{T- ^l^Tx^-K-) 


1.2 X 10"^ (e") 


2.0 X 10-'' (e+) 


2.9 X 10"'^ (/i") 


7 X 10-** (^+) 





Let us survey some impacts of the Qv^P constraints. Presently the best experimental lower bound on the Ot^/3/3-decay 
half life is that obtained for "^^Ge 0: 



T°/2(™Ge) > 1.9 X lO^Syrs. 



(44) 



In Ref. [3] this bound was used to constrain the contribution of Majorana neutrinos of arbitrary mass. In a good 
approximation this constrain reads: 



E 






< 5x 10-**GeV-^ 



(45) 



with go = 105 MeV. Here ak denotes the Majorana CP-phase of the Majorana neutrino state v^ of mass m^k- For 
the light-heavy neutrino scenario one has: 



E 



\u, 



eN\ 



Af=heavy 



rriN 



%^ E 1^- 

i=light 



< 5 X 10-^GeV"^ 



(46) 



where rriM ^ go ^-nd rriui <^ go- This upper limit, applied to each term separately, leads to a very stringent constraint 
on the heavy neutrino mass and mixing 



rriN 



<5x 10-^GeV-^ 



(47) 



One may use this limit to constraint LNV and LFV processes with e^ in the final state. From Eq. (32), with the fact 
that ai{mN)\UiN\'^ < Fjv, we obtain a useful inequality 



X- f^ \ - \UhN\'^\Ul2N\'^ . \UhN\'^ , , 

■FhhK^N) = t; < — : li,l2 = e,fj,,T 



N 



a-u 



valid for all masses mjv. Gombining Eqs. mm and (47) we get 



(m^lMeV) 11 

Fei < • 5 X 10 for I = e, /i, r. 

ai 



(48) 



(49) 
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FIG. 2: Exclusion plots for the sterile neutrino mass niN and mixing matrix element It/eivp from various experimental searches 
(for a summary and references see Ref. [IB]). The thick solid line is our exclusion curve, derived from the existing experimental 
upper bounds Ref. [16] on the rates of r and K, D, _B-meson semileptonic LNV and LFV decays, using the Monte Carlo 
sampling method. The inclusive approach is applied for the calculation of the total decay rate of a heavy Majorana neutrino 

iV. 



Using these limits in Eqs. (15)-(17l we can derive upper limits on some LNV and LFV decays imposed by non- 
observations of 0!//3/3-decay ( 44 1 . Varying itin within the resonant regions ( 14 ) of the considered decays we determine 



the largest rates of these decays, which we present in Table II. These limits represent absolute upper bounds on 
given decays within the sterile neutrino extension of the SM, compatible with Oi'/3/3-constraints. Here an important 
assumption is implied. We assumed that CP is conserved in the neutrino sector. This case corresponds to no 
cancellation between different terms in Eq. (45), when the limit (47) is valid. If there are more than one heavy 



neutrino state Ni~, with different CP violating phases ajv, then different terms in the first sum of (46 1 may compensate 
each other, reducing the individual Oi^/3/3-constrain on each of them. Note that even with the presence of only one 
heavy neutrino there may happen this kind of compensation between the heavy N and the three light neutrinos t'i, 
since they all coherently contribute to 0i//3f3-deca.y. Thus we conclude that observation of some of the processes in 
Table II, at rates above the indicated limits from Ov/ijS, shown in square brackets, may point to the following; (a) 
Presence of a sterile Majorana neutrino N with the mass ttin in the resonant range of these processes and that CP in 
the neutrino sector is violated; (b) Presence of a Dirac sterile neutrino in this mass range, if LFV decay is observed; 
(c) There is an exotic low mass scale new physics with light particles within the same resonant range of masses, but 
other than a sterile neutrino. So far most of the existing experimental limits shown in Table II are significantly weaker 
than the corresponding Oz^/3/3 limits, except for K — >■ cctt, iien. We note, however, that the latter is not yet definite 
and requires more detailed analysis 113] of the corresponding experimental data. Uncertainties are related to the fact 
that for small mixing only part of the produced heavy sterile neutrinos decay in the detector and as a consequence 
the actual experimental bounds on K ^ eevr, fxen could be several orders of magnitude lower than indicated in Table 
II. Therefore more careful analysis of the existing experimental data and further searching for decays with e* in 
final states may provide an important insight into CP violation in the neutrino sector. The experimental bounds on 
branching ratios of the considered decays are expected to improve in the near future by an order of magnitude or 
more. This is, in particular, possible in LHCb experiment. 
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FIG. 3: Exclusion plots for the sterile neutrino mass tun and mixing matrix element |[/p]vp from various experimental searches 
(for a summary and references see Ref. [I3]). The thick solid line is our exclusion curve derived from the existing experimental 
upper bounds Ref. [16] on the rates of t and K, D, _B-meson semileptonic LNV and LFV decays, using the Monte Carlo 
sampling method. The inclusive approach is applied for the calculation of the total decay rate of a heavy Majorana neutrino 

N. 



B. Limits on as yet experimentally unconstrained decays. 



Presently there are no experimental limits on semileptonic decays of mesons with t in final states Mi — > t ?2 M2, 
and on any semileptonic decay of B^ like B^ -^ h h M. Below we derive theoretical upper limits for these decays 
within the considered sterile neutrino extension of the SM, on the basis of the fact that the maximal rates are reached 
in the resonant range for the sterile neutrino masses. Although such limits are linked to a specific scenario, they may 
have a more general meaning, considering the fact that the sterile neutrino extension of the SM is the only reasonable 
model having a light particle, sterile neutrino iV, which can cause resonant enhancement of the processes in question. 
Other known beyond the SM scenarios are associated with larger mass scales and heavier particles outside the resonant 
regions of these processes and then their contribution is expected to be much smaller than that from the resonant 
mass sterile neutrino. 



For derivation of the limits on as yet experimentally unconstrained decays we use the constraints ( 49 1 and 



J'rl < 



r^l < 



7rG'(Zr)"T-Af 



J'lu < 



\U, 



TTKy (G'^(zm) + G^li{zM))mN 



(50) 



l-iN\ 

ai 



for I = fi,T] 



\UrN\^ < 0.018 for rriN > nir- 



The first two constraints are derived from Eqs. (15)-(17l and are valid in the range tti^ + mi < niN < rrir — to^ 
for r-decays and for meson decays in the ranges indicated in Eqs. (14). The third constraint is a particular case 



i^ij 



of the inequality (48). The last one originates from precision electro- weak measurements [SI I20j . In order to derive 
upper limits on as yet unconstrained decays we scan their corresponding resonant ranges (14 1, selecting for each value 
of rriN the largest value of their rates compatible with the constraints in Eqs. (49) and (50 1. In Eqs. (50 1 we use 
the experimental upper bounds T'^^p{t — >■ I tt tt), T'^^p{M 
select the most stringent of the known (see Fig. 



tt) shown in Table II. For |C/^Ar| in Eq. (50) we 



3) limits for a given value of ttijv. The resulting upper limits on 
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FIG. 4: Exclusion curves for the sterile neutrino mass mjv and mixing matrix elements \UenUtn\ and \U^nUtn\, derived from 
the existing experimental upper bounds Ref. [16] on the rates of r and K, D, B-meson semileptonic LNV and LFV decays using 
the Monte Carlo sampling method. The inclusive approach is applied for the calculation of the total decay rate of a heavy 
Majorana neutrino N . 

the branching ratios for some LFV and LNV processes are listed in Table II in curly brackets. Note that the limits 
presented for the decays with e* in the final states imply CP conservation in the neutrino sector, discussed in the 
previous subsection. As seen from Table II only the experiments searching for the LNV and LFV decays of D and 
Dg are rather close to probe the rates of these processes not excluded by our analysis. Searches for other decays are 
far away from this perspective. As to the iiT-meson LFV and LNV decays, we refer to our comments at the end of 
the previous subsection. 



VI. SUMMARY AND CONCLUSIONS 



We have studied LFV and LNV decays of r and B,D,K mesons mediated by heavy Majorana neutrinos. We 
focussed on the dominant mechanism via resonant Majorana neutrino production. 

On the basis of joint analysis of experimental limits on certain sets of LNV and LFV decays with inter- 
secting resonant regions we extracted upper bounds on the heavy Majorana neutrino mass m^r and mixing 
\UeN\AUfj.N\AUeNUrN\AUfj.NUrN\ without additional assumptions on these parameters. Our limits derived from 
the above mentioned decays are shown in Figs. [2]|4] As seen they are less stringent than the corresponding limits 
derived in the literature. To our best knowledge the latter typically involve various ad hoc assumptions on the relative 
size of the heavy neutrino mixing matrix elements. Therefore, it is natural that our limits, not referring to any kind 
of such assumptions, are less stringent. The mixing Utn is left unconstrained by the present experimental data, and 
requires searches for processes with tt in the final state. 

In our analysis we used the method of calculation of the total decay rate of heavy Majorana neutrino on the basis 
of an inclusive approach. Instead of a detailed calculation of the sum of all the possible decay channels with different 
hadrons in the final state, we approximated the hadronic final states by a gg-pair, as suggested by Bloom-Gilman 
duality. This allows us to avoid large theoretical uncertainties in the decay constants of heavy mesons. 
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Special emphasis has been made on the stringent hmits on the heavy Majorana neutrino from non-observation of 
Oi//3/3-decay, from which we also derived indirect limits on various LNV and LFV processes. Our results, displayed 
in Fig. [2] and Table II, have been obtained assuming no CP-violation in the neutrino sector. As seen from Table 
II, our expectations for the prospects of observation of the decays with both e^e^ and e^e^ in the final states are 
pessimistic. One of the messages of the present paper is that despite of such a discouraging prediction these processes 
are worth searching. The point is that any observation of such decays above the limits set by the CP conserving Oi'(3(3 
decay constraints would most likely point to the existence of heavy Majorana neutrinos in the resonant regions of 
these processes, and CP-violation in the neutrino sector, as explained in sec. V. 

We applied our limits on heavy Majorana neutrino mass ruN and its mixing Uin for a prediction of upper limits on 
the rates of as yet experimentally unconstrained LFV and LNV decays of B and D mesons. These limits are shown 
in Table II, and indicate that only the experiments searching for the LNV and LFV decays of D and D^ are rather 
close to probe the rates of these processes not excluded by our analysis. Searches for other decays are far away from 
this perspective. 
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